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ABSTRACT 

Numerical solutions of the isothermal Lane-Emden equation are presented, corresponding to self-gravitating 
gaseous cores embedded within a finite density envelope of overall cylindrical symmetry. These structures may 
be members of a fragmentation hierarchy proceeding from sheets, to filaments, to elongated, prolate clumps. The 
embedded solutions are the first of their kind, and as such represent a significant improvement upon the isolated 
cloud paradigm used almost exclusively by previous authors. The properties of the equilibria are in reasonable 
agreement with observations of dense molecular cores in star-forming clouds, despite the fact that there is only 
one free parameter in the models. We show that this parameter may be identified with the critical wavelength for 
instability in the parent filament. The implications of further fragmentation and the possible influence of magnetic 
fields are briefly discussed. 

Subject headings: ISM: clouds — ISM: structure — hydrodynamics — instabilities — Stars: formation 



L INTRODUCTION 

The present observational picture of star-forming regions 
conveys cloud structures considerably more varied than any 
one theoretical scenario can explain. Nevertheless, there are 
certain points of correspondence, when one restricts attention 
to the molecular, self-gravitating component of the interstellar 
medium. On scales of tens to hundreds of parsecs, both atomic 
and molecular gas clouds appear shell-like and filamentary 
(Scalo 1985; Kulkarni & Heiles 1988). Shells may be formed 
by several distinct dynamical processes: cloud-cloud collisions 
(Smith 1980), compressional shock waves from supernovae or 
OB stars (McCray & Kafatos 1987), and large-scale shocks as- 
sociated with spiral density waves (Roberts 1969) are but a few 
possibilities. Filamentary structure persists down to scales of 
several parsecs, to the regime of individual molecular clouds 
[Loren 1989, Nozawaet al 1991 (pOph); Heyeret al 1987, On- 
ishi et al 1996 (Taurus); Bally et al 1987, Tatematsu et al 1993 
(Orion)]. Each filament typically contains several distinct sub- 
condensations in close proximity to each other, some of which 
harbor infrared continuum sources (Onishi et al 1998). In some 
cases, the embedded clumps are spaced quite regularly along 
the filament (Schneider & Elmegreen 1979; Dutrey et al 1991). 
Thus, several authors have speculated that the formation of star- 
forming clumps proceeds via a hierarchical fragmentation pro- 
cess, in which filaments are formed out of larger structures, and 
then clumps out of the filaments (Schneider & Elmegreen 1979; 
Gaida, Ungerechts, & Winnewisser 1984; Hanawa et al 1993; 
Fiege & Pudritz 2000a). 

Focusing now on the cores themselves, maps in dense tracers 
such as NH3 and CS display roughly elliptical intensity con- 
tours, with a mean apparent major- to-minor axis ratio of around 
2 (Jijina, Myers, & Adams 1999). Statistical arguments ap- 
plied to the distribution of measured axial ratios for several sur- 
veys have prompted some authors to conclude that the cores 



are more likely to be intrinsically prolate than oblate (David 
& Verschueren 1987; Myers et al 1991; Ryden 1996). Other 
statistical arguments imply that the observed elongation is un- 
likely to be a result of star formation or outflows (Myers et al 
1991). Coupled with the additional result that the majority of 
cores are near virial equilibrium (Jijina et al 1999), this argues 
against a wholly dynamical origin for prolateness, since the 
implied lifetimes are then so short that observations of cores 
without embedded stars would be exceedingly rare (the latter 
represent about one-half of the ammonia cores detected in the 
Taurus region, and a somewhat smaller fraction in Ophiuchus 
and Orion). Rather, it appears more likely that cores are at least 
quasi-equilibrium structures, and that their shapes therefore of- 
fer some clue to the forces responsible for their formation. 

If this reasoning is correct, then it leads to a formidable crisis 
in our current theoretical picture of cloud equilibria, which typ- 
ically envisions star forming cores as self-gravitating clumps 
bounded by a zero density, constant pressure — and so high 
temperature — medium (e.g., McKee et al 1993). That is, it 
is difficult to conceive of prolate, quasi-equilibrium, gaseous 
cores as purely isolated structures. Non-gravitational forces, 
such as rotation or magnetic fields, are not likely to aid in main- 
taining prolate equilibria, although their role in oblate structures 
has been made abundantly clear^ The absence in the literature 
of even a single, physically acceptable, gaseous, prolate equi- 
librium solution suggests that a fresh theoretical approach is 
necessary. We examine a scenario in this paper whereby core 
morphology is directly attributable to a fragmentation process. 
As we show in §|| and ^ the existence of prolate structures 
within cylindrical filaments may be understood in exactly the 
same way as the existence of the filaments themselves: as a 
result of the fragmentation of the parent cloud and nonlinear 
growth of the fragments. What has apparently escaped attention 
until now, and what we demonstrate explicitly, is that accessi- 
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-It is in fact possible to construct isolated, prolate magnetic clouds, but these possess magnetic field structures which are highly unusual and, in any case, lack 
direct observational justification. See Fiege & Pudritz (2000^)) and Curry & Stabler (2000) for examples. 
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ble, long-lived states exist wherein cores and their extensive fil- 
amentary envelopes occupy the same hydrostatic structure. We 
now proceed to outline a particular fragmentation hierarchy that 
might be responsible for such structures. 

2. NONLINEAR FRAGMENTATION OF AN ISOTHERMAL LAYER 

To motivate the new solutions presented in the following sec- 
tion, we begin by outlining a possible formation mechanism for 
filamentary structures in the molecular interstellar medium. By 
drawing attention to selected existing results in the literature, 
we hope to provide a more compelling argument for a partic- 
ular fragmentation hierarchy alluded to previously by several 
authors, notably Schneider & Elmegreen (1979) and Larson 
(1985). 

2.1. Analytic Solutions 

While not realistic in detail, a large, shell-like structure in the 
interstellar medium may be approximated on sufficiently small 
scales by a planar, self-gravitating layer, or "sheet" (Elmegreen 
& Elmegreen 1978). We restrict consideration throughout to an 
isothermal, self-gravitating gas, whose equilibrium is governed 
by the Lane-Emden equation. 
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where is the gravitational potential, a is the constant sound 
speed, and po.s is a constant. The equation has been written 
in cartesian coordinates (jCjjjZ). When t/j varies solely in the 
y-direction, i.e. = V'id(3'), and under the usual boundary con- 
ditions ■0id(3' = 0) = 0, ((i'!/'iD/i^3')v=o = 0, an exact solution is 
known (Spitzer 1942): 

Pm{y) = POj exp[-?/;iD(}')/a^] = po,., sech^(y/£o), (2) 

where the scale height £o is defined by £Q = a/ (InGpo.sy^^. The 
layer is unbounded in the y-direction; solutions truncated by an 
external pressure at constant y may also be constructed (e.g., 
Elmegreen & Elmegreen 1978). However, as these introduce 
an additional parameter — namely, the external pressure — into 
the problem, we do not consider them here. The stability of 
the unbounded solution was examined by Ledoux (1951), who 
found gravitational instability for infinitesimal sinusoidal per- 
turbations of wavelength exceeding A i cr = 2tt£o. The maximum 
growth rate occurs at A.v.mgr = 2.24Av,ci- = 14.1^0. The choice of 
the ;ic-direction for the perturbation wavevector here is arbitrary; 
the same results hold for perturbations in the z-direction only. 

Subsequently, Schmid-Burgk (1967) (hereafter S-B) pre- 
sented a remarkable two-dimensional (2D) solution of equation 
(|l]) (independently discovered, in a hydrodynamic context, by 
Stuart 1967): 



P2D{x,y) = po.s exp 



-tp2D(x,y) 



Po.(l-A^) 



[cosh(3;/£o)-Acos (x/£q)]^' 
(3) 

The constant A, < A < 1 describes the amplitude of spatially 
periodic variations along one direction parallel to the layer (here 
chosen as x), with the same critical wavelength, X^.a, as found 
by Ledoux. As noted by S-B, this result shows that the latter 



mode is not restricted to infinitesimal amplitudes. Indeed, so- 
lutions of all amplitudes A are contained in equation ^? It 
is noteworthy that the periodicity of this exact, nonlinear solu- 
tion does not correspond to that of the fastest-growing mode in 
the linearly unstable layer, but rather to the critically unstable 
mode. This feature is discussed further in §^. 

As A ^ 0, equation (^ reduces to Spitzer's solution. In the 
limit A 1, the equilibrium'* corresponds to a series of parallel 
cylindrical "fragments," infinitely extended in the z-direction, 
with density maxima separated by a distance A^.a-. For general 
A, the fragments are embedded, nearly elliptic cylinders with 
equidensity contours of eccentricity (1-A)'/^ (see Figure 1 of 
S-B). The density maxima and minima along the x-axis (each 
spaced at an interval of Ax = lirio), are given by 



1+A 



and 



Ps = P{),s 



1-A 
T+A 



(4) 



respectively. 



2.2. Numerical Solution 



There is reason to suspect that analogues of the S-B solution 
exist in other geometries, in particular in cylindrical symmetry, 
a case we shall focus on in the following section. Unfortunately, 
to our knowledge, no other embedded solutions are known in 
analytic form. However, numerical techniques can be used. As 
a test of one such scheme used later in this paper, we attempted 
to find the S-B solution using the self-consistent field method 
(e.g. Tassoul 1978). That is, equation ([I|) was solved iteratively 
on a two-dimensional grid, subject to the boundary conditions 



dx' 

V''(0,0) = -lnp'(0,0), 



dj/_ 



and 



= 0, 

'=0 

a7 



(5) 



^0. (6) 



Here a prime indicates a dimensionless quantity — e.g., ip' = 
ijj/a'^ — and <x' <X', < / < F' is the extent of the com- 
putational region. The nondimensionalization of length will be 
discussed below; for now it may be assumed arbitrary. 

The first of conditions (^ ensures that the x-component of 
the gravitational field vanishes on the y-axis, while the second 
imposes reflection symmetry about the x = plane. Conditions 
(H) are consistent with the Spitzer solution, although neither 
constrains the numerical solution to be identical to the former 
in any particular limit. In addition, these conditions allow for 
solutions with a periodic structure in x. In an attempt to find 
these, we took ■(/;Jq ^iy') = -In ,,(3;') (see eq. ^ as an initial 
guess, and added a perturbation of the form 

Sip'(x',y') = e cos(px'/X') sm(py'/Y'), 

where e is a small, constant amplitude and Q < p < In. For con- 



l(n) l„lAn+\) 



< 



vergence of the numerical code we required l — ip"^">/ip 
S = 0.005, where ijj'^"^ is the value of i/j' at the n* iteration. For 
e < 0.005, the method converged immediately to the Spitzer 
solution (||). For any larger e, however, a unique, 2D struc- 
ture resulted. The equilibria so obtained constitute a family of 

^Note also that, like other solutions of the Lane-Emden equation. S-B's solution is homology invariant (Chandrasekhar 1958); i.e. the potential -(/"J]-, = 
ip2D(Cx,Cy)-2 log C, C = constant, is also a solution of equation (|l]). 

^The term "equilibrium" is used here in its strictly limited sense: i.e., a particular equilibrium may be either stable or unstable to perturbations in a direction in 
which the solution is translation-invariant (in this case, the z-direction). 
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solutions in the single parameter X', such that the density con- 
trast Pc/ps increases monotonically with increasing X'. Here 
Ps = p{X,Q). The vertical extent Y' was chosen sufficiently 
large that the effect of different Y' on 2D solutions of the same 
X' was negligible. 

In order to compare these solutions to the analytic one of S- 
B, we need to specify the nondimensionalization. This is, in 
fact, already implicit in the solution method. At each iteration, 
the Lane-Emden equation is solved for tp'ix' ,y') with the source 
term 

p'^"\x',y') = exp[ij'^"HO,0)-^'^"\x',y')]. 

Thus, p'{x',y') is normalized with respect to its central value at 
each iteration; i.e. /3'*"*(0,0) = 1 (this is made explicit in the first 
of equations ^). Consequently, once an exact, 2D equilibrium 
is found, all "memory" of the original density scale po j is lost. 
The characteristic density of the solution is pc = p(0,0), and 
thus the corresponding unit length is Ic = a/ilTrGpc)^/^, not 

-1/2 

£{) oc Pq / as in the analytic solution. Hence, the appropriate 
nondimensionalization is specified a posteriori as 



{x',y',z!) = {x/lc,y/lc,z/lc). 



(7) 



An exact relation between the parameter X' in the numerical 
problem, and the amplitude A in the analytic solution, may now 
be derived.^ Solving equation (||) for ip2D{x,y) and substituting 
the result into the second of equations (ra) gives 



4 \Ps 



1/4 



1+A 



1/2 



(8) 



The parameter X' represents one-half of the spacing between 
the embedded fragments; i.e., A^. = 2X' . In the analytic so- 
lution, the spacing of the fragments is constant, i.e. Ax/£o = 
Ax,cr/^o = 27r. There the fragment spacing stays the same while 
the amplitude A varies. Just as analytic solutions exist for all A, 
numerical solutions exist for all values of X' > tt. 

Figure 1 of S-B (1967) shows a solution with A = 0.17, 
whose corresponding density contrast is pd Ps = 1-987. In- 
serting this A into equation (^, one finds X' = 3.730. Using 
this X' as input for the numerical method then yields a solu- 
tion with Pc/ps = 1.981, a difference of 0.3 percent from S-B's 
analytic result. This solution is displayed in Figure 1 . The dis- 
crepancy between the analytic and numerical solutions remains 
small as X' increases above 3.73, and is still less than 1 per- 
cent for X' = 22, Pc/ Ps ~ 2400. At smaller X', the numerical 
and analytic results slowly diverge, with the relative error in 
Pc/ps reaching 5 percent at X' = 3.32, and nearly 20 percent 
at X' = TT. At the latter value, equation (^ predicts A = and 
Pe/ps = I- Presumably, the method is not sufficiently sensitive 
to the slight density contrasts present in these small amplitude 
(A < 0.05) solutions to render an accurate result. On the other 
hand, for X' < 3.125, the numerical method converges to the 
Spitzer solution, with a slightly non-uniform density along the 
midplane (e.g., Pc/ps — 0.99 at X' = 3.12). Overall, these re- 
sults confirm the validity of the numerical method, with the 
nondimensional length and density chosen as in equation (^. 
We may expect the same to be true for other solutions found by 
the same technique, iiTespective of whether an analytic solution 
is known. 



We conclude our discussion of the S-B solution by noting 
that related states have in fact appeared in the literature, in the 
context of either linear perturbation theory (Miyama, Narita, 
& Hayashi 1987fl) or time-dependent nonlinear calculations 
(Miyama, Narita, & Hayashi 1987^7). The basic character of 
the solution persists even in the presence of a (z-independent) 
magnetic field (Fleischer 1998; Nagai, Inutsuka, & Miyama 
1998), or when the layer is truncated by an external medium. 
The pressure-bounded, magnetized layer has additional modes 
of fragmentation available to it, as shown by Nagai et al (1998). 
However, the significance of the basic result — namely, that the 
embedded filaments are in fact exact, equilibrium solutions — is 
rarely emphasized. 

3. NONLINEAR FRAGMENTATION OF AN ISOTHERMAL FILAMENT 

3.1. Equilibrium and Stability 

In cylindrical coordinates (r,(j),z), the isothermal Lane- 
Emden equation reads 



(9) 



where poj is a constant. Here the z-axis is taken to co- 
incide with the principal axis of the cylinder, and we have 
taken d/dcj) = 0, restricting consideration to axisymmetric so- 
lutions. When ip depends only on r, i.e. = "ipiuir), and un- 
der the boundary conditions ■0id('" = 0) = 0, (diljiu/dr)r=() = 0, 
StodoUciewicz (1963) and Ostriker (1964) derived the following 
exact solution of equation (H) (hereafter the S-O solution): 



P\D(r) = Poj exp 



POJ 



(10) 



where Zq = {2a^ /irGpoj)^/^ is the cylindrical scale radius. In- 
terestingly, solution (|iy) reduces to that of S-B as A ^ 1 in the 
latter; details may be found in the Appendix. This suggests that 
equilibria quite similar to the S-O solution may result directly 
from the fragmentation of an isothermal layer. 

The solution ( p^ decreases as r~^ at large r; as in the pla- 
nar case, solutions truncated by an external pressure at constant 
r have also been considered (e.g., Inutsuka & Miyama 1997; 
Fiege & Pudritz 2000fl). The maximum mass per unit length of 
the filament described by equation (|^) is^ 



Pn 



2fl2 

2T:piu(r)rdr= —. 
^ 



(11) 



Isothermal cylinders with density distributions having p, > p^nax 
are unstable to radial collapse, while those with p < p^^^^ ex- 
pand radially outward, unless confined by an external pres- 
sure (e.g ., Inutsuka & Miyama 1992). The stability of solu- 
tion (ntt) to axisymmetric, linear perturbations was examined 
by Stodolkiewicz (1963), who found instability for perturba- 
tions of wavelength exceeding Aj ^i- = 3.94 /o- The maximum 
growth rate occurs at Aj.mgr = 1-98 A^-.d- = 7.82 Iq (Nagasawa 
1987). Further analysis has been carried out by Inutsuka & 
Miyama (1997) in the nonlinear regime of perturbation growth. 
However, the nonlinear resolution of the instability remains an 
open question. That is, what is the final outcome of the frag- 
mentation instability in an isothermal cylinder? 



^I thank Dean McLaughlin for pointing out a simplified derivation. 

*In the S-B solution, the total mass per unit z-length of each "cell", —n <x/Iq < n, is also 2a^/G. 
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3.2. Numerical Solution 

We undertake a numerical investigation of the question p osed 
above, using an identical technique to that described in §2.2. 
That is, rather than following the time evolution of a particu- 
lar unstable state, we instead search for an exact, static solution 
having a 2D structure. The boundary conditions on the poten- 
tial are now 



V''(0,0) = -lnp'(0,0), 



and 



= 0, 

'=() 

a? 



(12) 
^0,(13) 



where < r' < R', < z' < Z' is the size of the computa- 
tional region. The nondimensionalization is the same as that 
used in §2.2, except with 4- replaced by 1^ = (la^/irGpcy^^, 
where pc = p{Q,Q). Again we took the equilibrium potential, 
V''id('"') - p'id('"') initial guess, and added a perturba- 
tion of the form 

6ij'(r',z') = e sinipr'/R') cos(pz'/Z'), 

with Q < p < 2tt. For any e > 0.005, a family of unique 2D 
structures was again found, now parameterized by Z'. The cor- 
responding R' was chosen sufficiently large that its effect on the 
solutions was negligible. In practice, choosing R' to be approx- 
imately twice the radius of the "tidal lobe" — i.e., the last closed 
isodensity contour, of density = p{Q,Z') — was sufficient. 

The density contrast between the center and the tidal lobe, 
Pc/ps, is a monotonically increasing function of Z'. Mod- 
els with 1.2 < Pc/ps <^ were generated, corresponding to 
2.1 < Z' < 30. A few representative equilibria are displayed 
in Figure^. At small Z', the equidensity contours are highly 
prolate (Fig. Ija); as Z' increases, the interior contours become 
nearly spherical (Fig. ^). Note that we were able to find 2D 
equilibria only down to a minimum Z' ~ 2.1; below this value, 
convergence was not obtained until Z' = 1.8. For all Z' < 1.8, 
the method converged to the S-O solution, with a slightly non- 
uniform density along the axis of symmetry (e.g., pc/ Ps = 0.997 
at Z' = 1.8). As we discuss below in §^ these results strongly 
suggest a fragmentation scale ~ Acr, as in the planar layer 

Physically, the quantity 2Z' may be interpreted as the separa- 
tion between any two clumps in a linear chain of identical con- 
densations (e.g. Lizano & Shu 1989; Fiedler & Mouschovias 
1992). This feature thus mimics an observed property of star- 
forming environments: namely, that pre-stellar cores are rarely 
found in isolation. In §i.2, we examine how the properties of 
the condensations depend upon Z', and how the latter may be 
constrained by observations. 

4. PHYSICAL PROPERTIES OF THE FRAGMENTS 

4. 1 . Virial Theorem Analysis 

Insight into the global properties of equilibria may be ob- 
tained from the scalar virial theorem. In the present context of 
embedded cores, however, care must be taken regarding the re- 
gion of application. For an axisymmetric cloud in equilibrium, 
we have (McKee et al 1993) 



where 



2U + U + W = 0, 



3 f 3 
U=- PdV = -a^M, 



(14) 



and 



Pr dS, 



p r • VtpdV 



are the relevant thermal, compressive, and gravitational ener- 
gies. These integral quantities are usually summed over the 
volume V enclosed within an arbitrary closed surface S. In the 
present case, one might choose for S any closed isobar. Then 
the thermal and compressive terms are easily calculated. How- 
ever, for any choice of S the gravitational energy W is incom- 
plete because part of the gravitational field arises from outside 
the region considered. This is in keeping with the nature of 
embedded equilibria. Indeed, each equilibrium found by the 
numerical method of the previous section comprises an entire 
"cell," i.e., the region 0<r<R,Q<(j)< 2tt, -Z < z < Z, so 
that the above choice of S is inappropriate in any case. The 
virial theorem may only be consistently applied on the same 
region. Since the gravitational force vanishes on the plane cir- 
cular section z = ±Z (by the second of the boundary conditions 
eq. |l3|), the corresponding energy W summed over the entire 
cell between z = -Z and z = +Z is therefore complete. 

Figure 51 shows the behavior of individual terms in the virial 
equation ([[4[), expressed in dimensionless units, along with 
their sum, S = 2U + II + W, as one proceeds along a sequence 
of increasing central density Pc/ps- Here, is used as a con- 
venient reference density only; the virial terms have been cal- 
culated for the entire cell, not just for the region inside the tidal 
lobe. The figure shows that the compressive term 11 has a mag- 
nitude of only 40 to 50 percent of the gravitational term W over 
most of the sequence. This may be compared with the case 
of isolated clouds, where external compression exceeds the ef- 
fect of self-gravity at low density contrast, the latter becoming 
dominant only for more centrally-condensed clouds. There, the 
boundary pressure is all that prevents the cloud from dispersing. 
That there are no such equilibria in the present sequence can 
also be seen by considering the total energy, E = W + U, also 
plotted in Fig. 0. Since E <0 for all pc/ Ps, all of the equilibria 
are gravitationally dominated. Indeed, this feature is consistent 
with the notion that such structures are the result of nonlinear 
fragmentation within the self-gravitating parent filament. Note 
that this also means that the prolate shape of the cores is no? pri- 
marily due to tidal stretching by adjacent cores; an effect that, 
in the isolated case, may be mimicked by use of the boundary 
condition eq. ( |l3| ) (see, e.g.. Curry & Stabler 2000). Finally, 
note that all three of f/, 11 and W approach nonzero limiting 
values as pd Ps ~^ 1- This is to be expected, since these quan- 
tities are nonzero in the S-O cylinder Their limiting values are 
readily calculated as: f/o = 37rZ', Ho = -2ttZ' , and Wo = -47rZ'. 

4.2. A Sequence of Prolate Equilibria 

As indicated by the dashed curves in Fig. ^ for a given Z' 
there exists a unique tidal lobe, within which the equidensity 
countours are closed, and that corresponds to a density mini- 
mum along the symmetry axis, r = 0. The intersection of the 
tidal lobe and the z = plane is a circle whose radius we de- 
note by R,. The tidal lobe is thus one possible definition of the 
"surface" of the dense core, since it marks where the latter may 
be distinguished from the background filament. However, these 
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solutions make clear the danger of taking such a term too liter- 
ally. On this interpretation, we may fix the density on the tidal 
lobe at a value appropriate to the intercore medium. (Defining 
the extent of a core using a single equidensity contour is in fact 
in accord with the usual observational definition, which assigns 
a core's size on the basis of its associated half-power intensity 
contour). 

Fixing the intercore density ps = mn, and temperature T al- 
lows one to characterize the size and mass of the fragments via 
the following reference quantities: 



Is 



a 

Gps 



1/2 



2 Ps 



1/2 



Ir 



{G^PsY'^' 
(15) 

In terms of these reference values, we now define the dimen- 
sionless lengths and mass: 



z=- = z' 

Is \2 Ps 



-1/2 



~ R, ~ M, 

R, = -r, and M, = —, (16) 



where M, is the dimensional mass contained within the tidal 
lobe. 

Figure ^ shows the behavior of Z. and M, as a function of 
density contrast for fixed a and Wj. The results are reminiscent 
of the Bonnor-Ebert (hereafter B-E) sequence, shown by dashed 
lines, but with a maximum in mass at pc/ Ps — 10.0 instead of at 
(Pc/Ps)be = 14.04. The maximum mass of the prolate sequence, 
Mt = 1.32, exceeds that of the B-E sequence, Mbe = 1182 by 
12 percent. The polar radius Z decreases from a value of 1.53 
at Pc/ Pi ~ 1 .2 to a minimum of 0.60 at p^/ Ps — 600. The tidal 
lobe radius R, has a maximum of 0.45 at Pc/ps — 6.4, and a 
minimum of 0.34 at Pc/ps — lO"'. 

Finally, it is instructive to compare the density structure of 
the prolate fragments with that of known equilibria. The equa- 
torial and polar density profiles for the models displayed in 
Figs. ^ and d are shown in Figures [sja and b. These figures 
also show the radial density profiles of a S-O cylinder and of 
a marginally stable B-E sphere. At the relatively low Pc/ps of 
Fig. gp, the polar density profile is shallower than z~^, while 
the equatorial density resembles that of a B-E sphere out to 
about three times its "core radius," after which it steepens to 
the of the S-O profile. At high p,/p, (Fig. |fo), the similar- 
ity of both density profiles to that of the B-E sphere is evident. 
Thus, except at very low central concentrations or at distances 
far outside the tidal lobe, the core's density structure bears little 
resemblance to that of its parent S-O cylinder This is consis- 
tent with observations of globular filaments, which rarely imply 
density profiles steeper than (Alves et al 1998; Johnstone & 
Bally 1999; Lada, Alves, & Lada 1999). 

4.3. Remarks on Stability 

In the case of isolated, pressure-bounded clouds, the pres- 
ence of a maximum in cloud mass as one proceeds from smaller 
to larger Pc/ps signifies a transition from stable to unstable 
equilibria. However, this technique (often referred to as the 
"static method"; Tassoul 1978), cannot be carried over to the 
embedded cores characterized by the mass M,, since they are 
not complete equilibria (§ [4.l|^ Thus, the maximum in the M, 
vs. Pc/ps relation seen in Fig.^ while suggestive, has no direct 
bearing on the issue of stability. 

Another method used to investigate the stability of isolated 
clouds is that of the Gibbs free energy (Stabler 1983; Tomisaka 



et al 1988). Unfortunately, this method is also unlikely to suc- 
ceed here, for the following reason. Along a sequence of equi- 
libria of fixed temperature and surface pressure, the Gibbs free 
energy is a minimum at the critical stability point. However, 
in the present context, an equilibrium consists of an entire cell, 
which lacks a single, isobaric bounding surface. It is there- 
fore impossible to construct a sequence of cells of fixed a and 
ns, for which quantities such as mass have definite maxima. 
Consequently, one finds that the Gibbs free energy is a mono- 
tonically decreasing function of Pc/p.?; i.e., no minimum value 
is attained along the sequence. Thus, the free energy approach 
cannot be applied in the usual manner. Explicit, time-dependent 
calculations of the evolution of small perturbations within the 
equilibria may be required to illuminate this important issue. 

5. COMPARISON WITH OBSERVATIONS 

To compare with observations, it is useful to write Is and 



as 



4 = 0.382 pc 



m, = 3.19 M, 







0.19kms"V V103 cm-3 

3 



-1/2 



(17) 



0.19 km s" 



The scaling for «j is the minimum total number density, n = 
I. Inn-,, estimated from '^CO measurements of filamentary 
clouds (Nercessian et al 1988), while that for a is the ther- 
mal sound speed for T = 10 K and a mean molecular weight 
of OT = 2.33 niH- In a high mass star-forming region such as 
Orion, «, can be as high as ~ lO'* cm"-' (Dutrey et al 1993), 
while a oc r'/^ can exceed 0.40 km s"'. Nonthermal motions, 
which dominate thermal motions in clouds above scales ~ few 
X 0.1 pc, may be included in a schematic manner by replacing 
a with a constant velocity dispersion a. Representative values 
of a are given below. 

With the above scalings, the results of Fig. ^ correspond to 
the following ranges: 

Q<M/Mq <4.2, 
0.23 pc < Z < 0.58 pc, 
Q<R, <0.17pc. 

These ranges are in reasonable agreement with deduced values 
for cores observed in dense tracers, such as ammonia (Jijina et 
al 1999). This is true even though NH3 has a critical excitation 
density, ngx ^ 10"* cm"^, an order of magnitude larger than our 
assumed n.,. However, this effect is offset in equations ( [l7| ) and 
( p^ by the fact that some regions have T > 10 K and that all 
exhibit nonthermal motions. 

We now compare the intercore separations found in this sim- 
ple model with those observed in star-forming regions. Two 
values of Z may be singled out from Fig. Q as being of particu- 
lar interest. The first, denoted by Zpeak, corresponds to the peak 
in M, and has the value Zpeak = 0.94. The second notable value 
is the minimum in Z occuring near the right-hand side of Fig. ^ 
Zmin = 0.60. The dimensional value of each of these quantities 
depends upon both the intercore density and velocity dispersion 
through equations ( [T^ and (|T7|). Figure ^ shows both charac- 
teristic separations as a function of n, for various cr, where the 
latter is given in terms of aio = 0.19 km s~', the thermal sound 
speed in r = 10 K gas. The range of observed intercore sep- 
arations, 2Zobs, in each of three star-forming regions — Taurus, 



6 



Ophiuchus, and Orion — is indicated by vertical bars on the 
graph. These ranges were deduced from molecular line maps of 
cores embedded within filaments in each region. The character- 
istic intercore density was estimated as the observed density at 
a scale corresponding to the mean value of 2Zobs for that region. 
Further details on the data used may be found in the caption to 
Table 1. 

Figure ^ shows that the observed separations agree with 
2Zpeak for cr/fl 10 ~ 1-2.5 in Ophiuchus and Orion, and for 
cr/aio ~ 1.5-3.5 in Taurus. The corresponding ranges com- 
puted for 2Zinin are cr/flio ~ 1.5-4 in Ophiuchus and Orion, 
and a/a\Q ~ 2 — 5 in Taurus. These ranges may now be com- 
pared with the observed intercore velocity dispersion (Tobs in 
each region. In Taurus, the range computed according to 2Zpeak 
includes (Jobs, while in Ophiuchus and Orion, dobs lies above 
the range of theoretical values. On the other hand, the ranges 
computed for 2Zmin include (Tobs for all three regions. However, 
since the physical relevance of Z ^m is highly uncertain due to 
unresolved stability issues (§4.3), henceforth we focus on the 
results for Zpeak- 

The discrepancy obtained using Zpeak can be stated alterna- 
tively as follows: given the observed a and in both Ophi- 
uchus and Orion, the model appears to overestimate 2Zobs- 
These estimates, 2Zpred, are given in Table 1 and plotted as open 
circles in Fig. ^. However, this disagreement may have a rather 
simple explanation. The model-derived values assume that the 
parent filament is aligned exactly perpendicular to the line of 
sight — an unlikely circumstance. Generally, the filament will 
be inclined at some angle 90° - / to the line of sight, so that 
2Zpred is reduced to 2Zp,edCos /. To bring this into agreement 
with the mean 2Zobs then requires / = 50° in Ophiuchus and 
;■ = 53° in Orion. 

An independent measure of the intercore velocity dispersion 
may be obtained from the observed properties of filaments. Ac- 
cording to equation ( [Tl| ) with a in place of a, the maximum 
line mass of the parent filament from which the cores could 
have condensed is: 



Mmax =16.4 



fllO 



pc ' 



(19) 



Using available observations, we computed an approximate line 
mass, /iobs, for filaments in Taurus, Ophiuchus, and Orion, with 
the results given in Table 1. The corresponding Umax, obtained 
from equation ( [T^ assuming /iobs = Mmax, is also tabulated. In 
Taurus and Ophiuchus, a^ax and cjobs are in close agreement, 
indicating that the filaments themselves are near virial equilib- 
rium. In Orion, Cniax should not be compared with (Tobs, since 
the latter is obtained from C''^0 measurements, whereas /iobs is 
derived from '^CO observations (Bally et al 1987). Instead, we 
note that the corresponding linewidth, Ay('^CO) = 2.2 kms"', 
gives a dispersion of (Tobs('^CO) = 0.97 kms"', again in very 
good agreement with (T^ax- 

6. DISCUSSION 

If the effect of inclination alone does not account for the 
above discrepancy between 2Zpeak and 2Zpred, then this may in- 
dicate a more fundamental limitation of the model in its present 
form. The smaller observed intercore separation (for a given 
lis and (t) might be the result of further fragmentation of stable 
cores lying at the low-density end of the equilibrium sequence. 
Fig. ^ However, the resulting cores would have correspond- 
ingly smaller masses than those obtained above in §pl which 



are already in only marginal agreement with observations. 

Although this paper has not included magnetic effects, it is 
worth mentioning that, in the presence of a longitudinal mag- 
netic field B that decreases in strength outward from the parent 
cylinder axis, the critical wavelength for instability is signifi- 
cantly reduced (Nakamura, Hanawa, & Nakano 1993; Hanawa 
et al 1993). This occurs because the effective sound speed, 
Ceff = (a^- v^)'/^, where va = B/iA-nGpy^^ is the Alfven speed, 
is decreased from its non-magnetic value, thereby reducing 
thermal support against gravity in the longitudinal direction. 
If a suitable generalization of the solutions found here exists 
for magnetized, isothermal cylinders, then the constituent frag- 
ments should have a smaller spacing for a given density con- 
trast, perhaps offering more satisfactory agreement with obser- 
vations. 

What is the relation between the intercore spacing and the 
most unstable and critical wavelengths in the S-O cylinder? We 
saw in §|| that, in the case of the isothermal layer, it is the crit- 
ical wavelength that determines the inter-filament spacing in 
S-B's 2D solution. This feature merits further comment. Let 
us imagine the temporal growth of an unstable fragment, either 
filamentary or prolate, to proceed as follows. Fragmentation be- 
gins at the scale of Amgr, the wavelength of the fastest-growing 
linearly unstable mode of the parental gas distribution. As the 
fragment (defined as the material residing within the tidal lobe) 
grows into the nonlinear regime, it begins to attract, and be 
attracted by, surrounding fragments with the same properties. 
Thus, the entire chain of masses becomes more tightly bound 
as a result of the increased gravitational potential energy due to 
condensation, and the intercore spacing decreases. There is a 
limit to this process, however, at the scale of Ad , below which 
the parent cloud is stable even in the linear regime. 

The above conjecture is supported in the planar case by 
the fact that the gravitational potential energy per unit length 
per cell increases in absolute value with the amplitude A (S-B 
1967). In the cylinder, the same trend is observed in W vs. Pc/ps 
(Fig. ph. Unfortunately, the lack of a 2D analytic solution — and 
therefore an amplitude — in the cylindrical case precludes the 
derivation of an algebraic relation between the intercore sepa- 
ration and Aj;.cr- However, a close inspection of the numerical 
results suggests that it is again A,_cr that is of primary impor- 
tance in filament fragmentation. 



Recall from the results of §3.2 that we were able to find 2D 
equilibria only down to a minimum Z' ~ 2.1, where Pc/ps = 
1.2. Below this value, the method converged to the S-O solu- 
tion. Now, if the minimum length scale for fragmentation were 
set by the critical wavelength as in the planar case, then we 
would expect to find 2D solutions down to a minimum value 



of Z' = Z//o = A-c-/(2/o) ~ 1.97 (§3.1), with a corresponding 



Pc/ps — 1- Our minimum converged value of Z' = 2.1 is con- 
sistent with this expectation, given the intrinsic inaccuracy of 
the numerical method. Moreover, the fact that equilibria with 
Pc/Ps> 1-2 were found for 2.1 < Z' < 3.91 = Aj,_mgr/(2/o) cer- 
tainly argues against any possible significance of A,,mgr in the 
final equilibria. 

Should a subsequent stability analysis of these solutions re- 
veal that the mass peak in Fig. ^indeed signifies a stability tran- 
sition, then an interesting, although highly approximate, crite- 
rion for whether a given filament will form stars can be derived. 
The argument proceeds as follows. The mass peak of Fig. Hcor- 
responds to Zpeak = 0.94, whence Zpeak = 0.94 4 = 0.36 (a/aio) 
pc. Thus, 2Zpeak = 0.72 (cr/aio) pc represents the scale of a sin- 
gle whole fragment formed out of the parent cloud. But what 
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is the affected region of the original filament? It is likely to 
be somewhat larger than 2Zpeak, by the argument summarized 
in the paragraph above. Roughly speaking, the affected region 
should be larger by a factor ~ Kmgr/ K.a- — 2, since A, „ gives 
the scale at which the filament first becomes unstable to frag- 
mentation (§3.1 ). Thus, assuming that the original filament was 
not radially collapsing (/i < fim^x), we estimate that any fila- 
ment longer than ~ 2 ■ 2Zpeak — 1-44 (cr/a lo) pc may be ex- 
pected to form stars via fragmentation and subsequent dynami- 
cal collapse of the prolate fragments. Shorter filaments, which 
could only harbor cores with pc/ Ps < 10, would require another 
agent — e.g., an increase in external pressure, ambipolar diffu- 
sion, etc. — to initiate collapse. 

7. CONCLUSION 

We have presented a new family of 2D, numerical solutions 
of the isothermal Lane-Emden equation in cylindrical sym- 
metry. The equilibria have an embedded, periodic structure, 
and in this sense are direct counterparts to the 2D solution of 



Schmid-Burgk (1967) in planar symmetry. Together, the ex- 
istence of the two equilibria suggest a fragmentation scheme 
which seems consistent with the observed hierarchical struc- 
ture in several well known star-forming regions. Moreover, our 
results constitute a remarkably simple yet robust explanation 
for the origin and maintenance of the prolate, gaseous cores 
that represent the lower rung of this hierarchy. Although we 
have ignored many physical effects that should be included in 
subsequent studies (particularly nonthermal motions and mag- 
netic fields), we expect that the shift in emphasis from isolated 
to embedded structures will prove fruitful in future theoretical 
work on the origin and evolution of dense cores. 

It is a pleasure to thank Steve Stabler for discussions that 
prompted a search for these solutions. 1 am grateful to Richard 
Larson and Dean McLaughlin for useful comments on the 
manuscript, and to an anonymous referee for suggestions which 
helped to clarify certain properties of the equilibria. 



APPENDIX 



RELATION BETWEEN THE SOLUTIONS OF S-B AND S-O 



It is worthwhile noting an intriguing correspondence between the solutions of S-B and S-O. Let us examine the structure of one of 
the cylindrical fragments of §^in its own right. That is, consider the limit A ^ 1, so that po,.? oc £q^ 0, giving an isolated cylinder 
(see equation [Q]). Let ^y' = 27r/Ax,ci- = ^.v.cr- Then in the limit fci cr the solution (^ becomes 



P(x,y) = Pc 



(Al) 



In terms of k-^.a, equation dS) reads 



the substitution of which into equation (|A1[) gives 



l-A 
T+A' 



p(x,y) = Pc 



= Pc 



Ax+3^ 4 
£2(1+A) ' 

^2 



2 el 



+ OiklJ 



(A2) 



in the limit as A 1, where we have noted that = x^+y^, appropriate to the geometry of the isolated cylindrical fragment. 
The S-O solution, on the other hand, reduces to the following in the limit of small r (equation [ITr 



p{r,z) = pQ.f 



1- 



2_r2 
ll 



'0 



(A3) 



Equations ( K2[ ) and (A3) are identical, at order r^, provided that: (i) /q = 4^^; and (ii) po.f = Pc- Recalling the definitions of and 
l[) given in §g and ^ if the latter is true, then so is the former. But equality (ii) does indeed hold, since the planar fragment has a 
topology indistinguishable from that of a cylinder in the prescribed limit. Thus, the S-O solution appears as the unique A ^ 1 limit 
of the S-B solution. 
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Fig. 1. — The closest numerical equivalent to S-B's A = 0.17 solution. Solid curves are equidensity contours, linearly spaced in p from p(Q,0)/ps 
p/ps =0.18. The dashed curve indicates the surface p = ps. Recall thatx' =x/£c,- ■ ■ etc. 



= 1.98 down to 
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Fig. 2.— Four embedded cores in the isothermal filament. Each is characterized by its value of (Z',R'): (a) (2.1, 2.1); (b) (3.0, 5.0); (c) (3.79,5.0); (d) (25.0,28.0). 
The corresponding density contrasts, Pc/ps, are 1.17, 4.70, 10.49, and 1076, respectively. Contours are spaced Unearly in (a), logarithmically in (b), (c), and (d). 
Recall that r' = r/fc, . . . etc. 
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Fig. 3. — Solid curves represent individual terms of the virial theorem for a n en tire cell, plotted against core density contrast (center to tidal lobe). The dotted 
curve is the sum of these terms, while the dashed curve is the total energy. See 34. 1 for details. 
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Fig. 4. — Solid curves show the entire sequence of models that were calculated, represented as Z.R,, and M, vs. central concentration, Pcj ps' The overall behavior 
is similar to that of the B-E sequence, shown by dashed hnes (dimensionless mass, top; dimensionless radius, bottom). See text for details. 
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Fig. 5. — (a) Density profiles along the polar (z') and equatorial (r') axes of the equilibrium depicted in (c)j|>f Fig. Q. Also shown are the radial profiles of a S-0 
cylinder (dotted curve) and a B-E sphere (dashed curve), (b) Same as (a), for the equilibrium depicted in Fig. 0(d). 
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3 3.5 4 4.5 5 

logio ^3 (cm-3) 

Fig. 6. — Intercore separation, 2Z, as a function of nj. Solid curves are 2Z^^\ from bottom to top, u/oio = 1, 2, 3 and 4, respectively. Dashed ciu*ves are 2Ztjun; 
from bottom to top, c/uiq = 1, 2, 3, 4 and 5, respectively. Vertical segments indicate approximate ranges of 2Z observed in three regions: from left to right, Taurus, 
Ophiuchus, and Orion. Open circles are predicted values of 2Z, assuming no inclination of the filament to the line of sight. See Table 1 and text for details. 
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Table 1: Mean Properties of Selected Filaments in Taurus, Ophiuchus, and Orion. 



Region, 
Cloud 


fJ-ohs 


lOglO 'Is,obs 

(cm-3) 


CT'obs 

(kms*') 


2Zobs 

(pc) 


(pc) 




Refs. 
(kms"') 


Taurus, 
L1495 


100 


3.48 


0.47 


0.6—1.4 


1.02 


0.47 


1,2; 2; 3; 4 


Ophiuchus, 
L204 


110 


3.65 


0.48 


0.36—0.73 


0.86 


0.49 


5; 5; 5; 6 


Orion, 
L1641 


385 


4.28 


0.58 


0.18—0.42 


0.50 


0.92 


7; 8; 8; 9 



Note. — The filaments in each region were chosen on the basis of their large apparent axial ratios (length, L/width ^ 20) as observed in '^CO maps, from which 

both M and L, and thus fi, were deduced. The intercore separations, 2Zobs, are nearest core distances derived from C'*0 maps in Taurus and Ophiuchus and NH3 maps 
in Orion. The mean intercore density, n, obs, is the density corresponding to the mean 2Z„^,^ in each region. This requires results from the appropriate (scale-dependent) 
molecular tracer: "CO in Taurus and Ophiuchus, and C'^O in Orion. The velocity dispersion, (T(,(,s, is chosen in the same manner References to the literature are 
ordered as follows. The first reference in each row is the source of ^obs, followed by references to the sources of n,.obs, fobs, and 2Zobs, each separated by a semicolon. 

References.— (1) Gaida et al 1984; (2) Mizuno et al 1995; (3) Nercessian et al 1988; (4) Onishi et al 1996; (5) Nozawa et al 1991; (6) Tachihara et al 2000; (7) 
Bally et al 1987; (8) Dutrey et al 1993; (9) Cesaroni & Wilson 1994. 



